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Abstract 

The theoretical model for the steady state Mach number of electrostatic shocks formed in the 
interaction of two plasma slabs of arbitrary density and temperature [l| is generalized for relativistic 
electron temperatures. We find that the relativistic correction leads to lower Mach numbers, and 
as a consequence, ions are reflected with lower energies. The steady state bulk velocity of the 
downstream population is introduced as an additional parameter to describe the transition between 
the minimum and maximum Mach numbers in dependence of the initial density and temperature 
ratios. Ion trapping and Landau damping are qualitatively discussed as dissipation mechanisms 
which can transform the soliton-like solution in the upstream into a shock. 
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I. INTRODUCTION 



Collisionless shocks can efficiently accelerate charged particles to high energies and its 
study is of interest to a wide range of scenarios, e. g. space and astrophysics, especially 
in the context of acceleration of cosmic rays to ultra high energies (up to ~ 10 21 eV) [2J, 
laser-plasma interactions with applications in proton therapy of tumors 31 or injection of 
particles for conventional accelerators 4j, and inertial confinement fusion |5j. Electrostatic 
shocks have been observed* interp.anetary space H, in parfcn.ar in the .onosphere Q 
and in the auroral region [10], where they arise from plasma cloud collisions at the interaction 
between solar wind and magnetosphere or between solar wind and interstellar medium in 
the heliosphere region 

Recently, electrostatic shocks have been propelled into the focus of research due to the 
ability of generating high Mach numbers in compact laboratory systems [3| , which could 
provide an alternative to the costly standard synchrotron accelerators. It was found that ions 
are efficiently accelerated in electrostatic shocks by reflection from the electrostatic potential 
with twice the shock velocity in the rest frame of the upstream [l^j]. Shocks with moderate 
Mach numbers were generated in a laboratory experiment, where ion beams with 20 MeV 
were produced 12]. The trend towards > 250 MeV/nucleon with a quasi-monoenergetic 
profile, which is relevant for the treatment of deep seated tumors 15] , was demonstrated by 
particle-in-cell simulations of laser-driven shock acceleration 13]. 

The theoretical framework for the interaction of two plasma slabs leading to non-linear 



structures with large Mach numbers [16|, [17] has been generalized by Sorasio et al. [l| for 



arbitrary plasma temperatures and densities. In this paper, we extend this work to relativis- 
tic electron temperatures and discuss efficient acceleration of ions.^he conditions for shock 
solutions are derived from an analysis of the Sagdeev potential 



being fully described 



by the ratios of the initial densities and temperatures, which also determine the condition 
for reflection and subsequent acceleration of ions by the electrostatic shock potential. 

In section [III the theoretical framework is described and the conditions for obtaining 
a shock-like solution are given. The dependence of the steady-state shock Mach number 
on the initial plasma temperature and density ratios is analyzed in section 1111} and the 
underlying characteristics leading to a minimum and a maximum value for the Mach number 
are investigated, considering the impact of a downstream velocity on the shock Mach number. 
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The conditions to obtain ion reflection and high energy ions are addressed in section HVl and 
dissipative effects that can create the oscillatory downstream potential are discussed in 
section [V] The results are summarized in section ED 

II. ELECTROSTATIC SHOCK SOLUTIONS FROM ARBITRARY UPSTREAM 
AND DOWNSTREAM DENSITIES AND TEMPERATURES 

Plasmas composed of hot electrons and cold ions are governed by ion sound waves with 
phase velocity oj/k ~ ^/k^Tjrrii = c s . The high ion inertia creates a restoring force to 
the thermally expanding electrons and wave steepening due to non-linear effects can lead 
to the generation of electrostatic shocks. Dissipative effects in the downstream region are 
necessary to have laminar shock solutions, in contrast to a solitary wave, and are addressed 
in section |Vj In this section, we focus on the physics in the upstream region, making use of 
the Sagdeev formalism [l^ in order to determine the electrostatic potential. 

The collision of two semi-infinite plasma slabs with arbitrary temperature and density, 
will lead to the formation of an electrostatic shock if the fluid velocity is small compared 
to the thermal velocity of the electrons. In the rest frame of the shock, the upstream 
population of electrons and ions is moving with v s h towards the shock. An electrostatic 
potential is formed due to the different inertia of the particles, which is steady in time in the 
shock rest frame, and increases monotonically from 0o at x = xq in the far upstream until 
it reaches its maximum fa at x — x±, where the transition to the oscillatory downstream 
region is defined (see Fig. [T|). The double layer is then maintained by four populations. 
Two populations of free electrons in the upstream and downstream regions moving towards 
the shock, with kinetic energies higher than the potential energy, a population of trapped 
electrons whose kinetic energy is less than the potential energy and a population of cold 
ions flowing from the upstream to the downstream region while being decelerated by the 
electrostatic potential. 

The electron populations are treated kinetically, while the ions are treated as a fluid. 
The electron distributions have to be solution of the stationary Vlasov equation and can be 
determined knowing the particle distributions in the unperturbed plasma. The free (left) 
electron population, propagating from the upstream to the downstream region, is described 
by the relativistic drifting ID Juttner distribution function /o(7o) — ^ r oK]~ 1 (/io) 7o(7o — 
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FIG. 1: Electrostatic shock formation from the interaction of free and trapped electrons and a 
population of cold ions. The electrostatic potential shows a monotonous increase from <f>o in the 
upstream region (x < xq) to cfii in the downstream (x > x\), where it becomes oscillatory. 



I) -1 / 2 exp[— /io7o(l — PoPsh)] 19| with the normalized velocity of the electrons f3 = v /c > 0, 
the Lorentz factor 7 = (1 — /3q) -1 / 2 , the thermal parameter /i = m e c 2 /IcbTq, where T 
is the electron temperature, m e the electron mass and ks the Boltzmann constant, and 
normalized fluid velocity f3 s h = v s h/c. The normalization constant contains the density of 
the left electron population N in the upstream region (x < Xq) and the modified Bessel 
function of the second kind K\. In the limit of nonrelativistic electron temperatures, /^o 3> 1, 
the distribution function is approximated by a Maxwell-Boltzmann distribution and the 
thermal velocity can be introduced as v t h,o = \Jk^o[m e = c/^/po- The free electrons 
in the downstream region (right population) are described by /i(7i) = A'" 1 K7 1 (/ii) 7i(7 2 — 
l) _1 / 2 exp — /Ui7i(l — fiifid) + fc^^ with /?i = v\jc < 0, where the parameters have the 
same meaning as for the left population and are indicated with a subscript 1, and (3d is the 



fluid velocity of the downstream. To faci 
we stick to the notation of Sorasio et al. 



itate comparison with the nonrelativistic model 
1] and have multiplied the distribution function 
by a factor containing the potential step <fti — 4>q. Ni represents the density of the right 
electron population in the far upstream and sums up with the density of the left population 
to the ion density Ni = No + N± at x < xo to guarantee charge neutrality. The trapped 
electrons are represented by the flat-top distribution function in the relativistic notation 
fit = A r 1 K 1 " 1 (/i 1 ) 7i(7 2 — I) -1 / 2 exp(— /ii), according to the so-called "maximum-density- 



trapping" approximation [2fJ,|2l|, which guarantees /i(7i = j c ) = fit at the critical Lorentz 



factor 7 C = 1 + e(0i — (po)/m e c that discriminates between free ((3\ < —/3 C ) and trapped 
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electrons < C ), shown in Fig. [2] where the more convenient parameter u = P'j has 
been introduced. Since the fluid velocities /3 s h and fid are small compared to the thermal 
velocities, we will neglect this dependence in the following calculations. 

We introduce the electron Lorentz factor 7 e which accounts for the electrostatic 
potential in the shock frame and make use of the conservation of energy to write 
the upstream and downstream Lorentz factors as 7 e = 70 + e(0 — <f)o)/m e c 2 = 
71 — e(<pi — (p)/m e c 2 > 1. The electron density can then be computed as 
n e = J-j 00 feipfe) d^e, obtaining the electron densities in the upstream region 
uq(Alp) = N K^ 1 (fj, )e Alf f™ A 1 e" Mo7e 7 e (7g — l)^ 1 / 2 d^ e and downstream region ni(A(^) = 



jv rKrW©) 



- M o7e/o 7e(T 2 _ y-1/2 dle + 2e-^/ v /(l + A^//i O ) 2 -l 



with the dimensionless quantities Atp = e(<f) — <fto)no/ m e c 2 , T = Ni/Nq and = ^0/^1 

Using the fluid equations for ion mass and energy conservation and assuming that the 
ions are cold and that none of them is reflected at the potential, the ion density can be 
determined as rij(Ay?) = iVj / a/1 — 2 Aip / M 2 , where M = v s h/c s is the ion acoustic Mach 
number, c s = ^kBTo/mi the upstream ion sound speed and is the ion mass. The particle 
densities are then combined in Poisson's equation 

d 2 Ay? 



rii(Atp) n (A(p) + ni(A^) 



N; 



(1) 



where the normalized quantity x — x/Xd with Xd = a/ KbT / (4ire 2 Ni) has been introduced. 
Defining the right hand side of equation ([TJ) as — dty (A(p) / dAtp , the similarity to the har- 



" U c U o 



FIG. 2: Electron distribution functions upstream (/q ~ dashed) and downstream (/x + fit ~ solid) 
for two different electron temperatures fio = 5 (blue) and 50 (black) with T = 3, = 2, e((pi — 
<f)o)/m e c 2 = 2 and /3 sh = 0.02. 



5 



monic oscillator can be immediately noticed and equation ([T]) can be seen as the motion of a 
pseudo-particle in the Sagdeev potential \&(A<£>) [l8j, allowing to identify bounded solutions 
as possible shock solutions. Integration of equation (TT]) with respect to Aip leads to 

1 fdAcpV 



2 [~&x) + ^ Av? ) = COnst =: ^° 
with \l/ = ^(v^O; M, T, 0) and the non-linear Sagdeev potential given by 

tt(Ap, M, T, 6) = Pi(A<p, M) - P e0 (Aip, r) - P el (A v , V, 6) 



(2) 



(3) 



where the quantities Pi, P e o and P e \ represent the ion, the upstream electron and the down- 
stream electron pressures, respectively, which are defined as 



Pi (Ay?, M) = M 2 1 



P e0 (Ay2, r,/x ) 

p el (A^,r,e,/i ) 



l + r 

re 
l + r 



2A<^ 

1^" 



Ki(fio) J i 



c?7e~ 



7 + —) -1-1 



(4) 
(5) 



o?7e" 



-moCt-i)/©, 



+ s V s 2 — 1 — log s + Vs 2 — 1 



with s = 1 + Aip/i^Q. In the case of nonrelativistic temperatures 



(6) 

u,q ^> 1, equations ([5]) and 



OH]) can be integrated analytically, retrieving the expressions of ll| 
P-(A^,r) = -l-^^ + e^erfcv^-l) 



(7) 



p e T(A^,r,e) 



( 2 /A^ ^ erfc /A* 8 /A^ _ 



with erfc the complimentary error function. Note that the explicit dependence on the up- 
stream and downstream temperatures vanishes in this approximation. For highly relativistic 
electron temperatures, /x «C 1, the pressures are approximated by 

Atp(l - fi ) 



Pe r !(A^r,e) 



l + r 

ApT 

e(i + r) 



6 1 



e 



(9) 
(10) 



and the explicit dependence on the temperatures is still maintained. Figure |3] shows the 
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FIG. 3: Sagdeev potential V(<p) obtained from equations ©-© for M = 1.7, T = 3, G = 1 
and no = 0.1 (red), 1 (green), 5 (blue), 100 (orange). Represented by dashed lines are the non- 
relativistic and highly relativistic approximations, given by equations ([?])- (TlOj) . 

Sagdeev potential, given by equations (J3])-([S]), for upstream electron temperatures /x = 
0.1 — 100 for a constant Mach number M and a comparison with the nonrelativistic and 
highly relativistic approximations (17])- fll0j) . A higher temperature leads to larger absolute 
values of the minimum of which results in lower values for the Mach number at which 
ion reflection sets in, as will be shown in the following section. 

III. MACH NUMBER DEPENDENCE ON INITIAL PARAMETERS 

As can be seen from equation (J4]), the model holds for Aip < M 2 /2 := Aip cr . The ion 
pressure becomes imaginary when the electrostatic potential exceeds the ion kinetic energy 

eA0 > ^rriivl (11) 

and the ions are reflected by the shock potential. We define the Mach number at which 
ion reflection sets in as the maximum Mach number M max . In order to determine possible 
shock solutions with M max , we use equation (T5]) which gives the condition for the existence 
of a monotonic double layer solution as \& := ty(Aip, M, T, 0) — < 0. For a given Mach 
number M, a soliton-like solution is possible only if the electron pressure exceeds the ion 
pressure. The solutions are found numerically by solving \l/(M 2 /2, M, T, O) = and are 
shown in figure HI 

As already found in jlj , the analytical dependence of the maximum Mach number in the 




FIG. 4: Maximum Mach number versus density ratio for the highly relativistic case /io = 0.1 
(black) given by Eq. (|12p and the non-relativistic case given in l| (blue) for a temperature ratio 
6 = 1. The dependences for very small and very large density ratios are indicated by the dashed 
lines. 

non-relativistic approximation is given by M max = 3y / 7r©/8(l + r)/T, which is M max « 
3-y/7r0/8 for large density ratios and has a M max oc T" 1 dependence for low density ratios. 
Here, we find for the case of highly relativistic temperatures, /iQ<l, 



M - = f e ( 1 + ?(T^e)> < 12 > 

which is displayed in figure H] together with the non-relativistic expression. It can be 
easily seen that the maximum Mach number is constant for high density ratios as in 
the non-relativistic case, M max « V2@ for Y ^> 1, and has a dependence M max « 
v/20(l + /i O )/r(l -/x /e) oc r- 1 / 2 for r < l. The comparison of the non-relativistic 
and highly relativistic cases in figure H] for a temperature ratio G = 1 shows that for higher 
upstream electron temperatures the maximum Mach number is reduced. 

We analyze now the lower limit and the range of possible Mach numbers for given tem- 
perature and density ratios. The shape of the Sagdeev potential and thus the existence of 
shock solutions depends on the choice of Y and 9 and we can distinguish three different 
types of solutions which are shown in Figure [5j Case (1) represents the case where shock 
solutions exist for ^ = \P — < and Aip > 0. While the monotonously growing Sagdeev 
potential in case (3) does not allow for shock solutions, case (2) defines the threshold with 
Acp = and provides the conditions to determine the minimum Mach number, which are 
given by d^/dAip = and ty(Aip) = 0. While in the highly relativistic limit M = 1 is the 
lower limit, in the non-relativistic case a lower limit M > 1 exists. The Sagdeev potential 
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FIG. 5: Different types of the Sagdeev potential *(A<p) = * - 
is expanded for Aip <C 1 since we are looking for solutions A<p — > 0, obtaining 

*(A^,M,r,e) « a^ 2 



2M 2 (1 - 




,(13) 



which is a function of the upstream potential <po- The minimum Mach number can then be 
found by solving ty(A(p, M(tpo), V, 0) = with the Mach number at the minimum of the 
Sagdeev potential given by 

r r /2 

M{<po) = JWol { 1 - T T 1 TT2 ) ■ (14) 

[ [e^erfc[ v ^] + r(^4 v ^ + e^/e er f c [ v ^7e]jj J 

For large temperature ratios we find small deviations from the approximation ipo = 
which is equivalent to ^/q = in equation fl2]), see Figure |6] (a). Panel (b) shows the respective 
Sagdeev potentials with = for the minimum Mach number according to Ac/? = in 
black and for the approximated model with ip = in red. The exact solution allows for the 
formation of electrostatic shocks at slightly lower Mach numbers. 

The transition between minimum and maximum Mach number can be described as a 
function of the steady state ion speed in the downstream region, Vd- In the rest frame of 
the shock, the upstream ions propagate towards the shock with velocity v itUS := v s h and are 
decelerated by the shock potential cp to velocities in the downstream < Vi^ s < v s h. The 
velocity is v^ds = if the ions are completely stopped by the potential and v^ds = v s h if 
they are unaffected and stream freely in the downstream region. In the upstream frame this 
corresponds to ion downstream velocities —Vi < v^du '■= —Vd < 0. Starting once more from 
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FIG. 6: (a) M m i n (dashed) and M max (solid) for the approximation (^o = (red) and comparison 
with the exact solution (black) for T = 1. The dashed vertical line shows the transition from M m j n 
to M max for downstream velocities < Vd < v s h- (b) Sagdeev potential for the minimum Mach 
number with Aip = (black) and comparison with the approximation ipo = (red). 

the energy conservation for ions, we can relate the downstream ion speed directly with the 
shock potential Vd/c s = M — M 2 — 2 Aip and use 

r 

Vsh, 



M 2 



1-11-^ 



(15) 



to find the zeros of the Sagdeev potential ^(Atpd, M, T, G) to determine the Mach number 
M. This transition is shown in figure EK- When the shock propagates with a speed slightly 
above the minimum Mach number, the downstream will have almost the same speed as 
the upstream population due to the small potential jump that has only a weak effect on 
the particles. At the maximum Mach number, the potential jump is so strong that the 
downstream propagates with the same speed as the shock front. 



IV. ION REFLECTION FROM THE ELECTROSTATIC POTENTIAL 

In the context of applications of shock-accelerated ions, e. g. for medical purposes, it 
is important to achieve high energies and the realization of the theory is restricted by 
experimental feasibility. As shown in Fiuza et al. (2012) Jj3], for typical experimental 
conditions associated with laser-driven shocks, ion acceleration occurs close to the critical 
Mach number M ~ M cr , so that the reflected ions will have a velocity Ui re /i = ^ v sh — 2M cr c s , 
which is determined by the condition for the critical Mach number. To simply reflect ions 
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from the shock, figures H] and [6^ show that large density ratios T and low temperature ratios 
are favorable, since low Mach number shocks are easier to drive. 

In order to achieve high energy ions, large Mach numbers M cr and/or large ion sound 
speeds c s are needed. An increase in the Mach number can be gained by a high initial 
temperature ratio 0, which is equivalent with increasing the energy of the slabs for a fixed 
upstream c s , and a high piston velocity, which provides a high momentum transfer to the 



plasma (i.e. initial relative fluid velocity between the two slabs) 22]. The sound speed can 
be increased by increasing the actual value of the electron temperature. This works well in 
near-critical density plasmas as a significant fraction of the laser energy can be absorbed by 



the plasma [13|. 



V. DISSIPATIVE EFFECTS LEADING TO AN OSCILLATORY DOWNSTREAM 
POTENTIAL 

So far, we have described the solitary solution in the upstream region and neglected the 
processes leading to a shock solution. The inclusion of a constant population of trapped 
electrons in the downstream region is not enough to break the symmetry. In this section, 
we qualitatively discuss possible mechanisms that lead to shock formation, which are ion 
trapping in the downstream of the shock and electron Landau damping. 

It is assumed that ions are reflected from the shock front, which leads to an electro- 
statically stable solution only if the fraction of reflected ions is small and wave-particle 
interactions are negligible. Moreover, ion trapping in the downstream region can also be 
important. The electrostatic potential in the upstream region is then given by the soliton 
solution derived in the previous section. The downstream potential is described by a second 
Sagdeev potential, which considers the physics in this region, by applying for example the 
Bernstein-Greene-Kruskal (BGK) method 



where the oscillatory electrostatic potential 



is separated into individual sections with each section considering the trapping of a different 
species. Figure [7] shows a sketch of the ion phase space (v,x) to illustrate this approach. 
Ions entering the shock region from the upstream are decelerated to v < Vi, and can en- 
ter the downstream region only if the initial kinetic energy exceeds the electrostatic energy. 
Trapped ions in the downstream have a kinetic energy smaller than the potential energy and 
their phase spaces are closed. To derive the Sagdeev potential, the downstream is split into 
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FIG. 7: Ion phase space showing free and trapped populations 



several sections and the energy distribution of trapped and untrapped ions and electrons is 
considered in each region. 



If the fraction of reflected ions is sma 
actual shape of the downstream potential 



L the upstream potential is not affected by the 



241 ] . Using the analogy of the harmonic oscillator, 



the electrostatic potential swings in the first Sagdeev potential until it reaches the maximum 
4>max, entering then the second Sagdeev potential, where it does the continuous oscillatory 
motion. The conditions to match the potentials require therefore \Pj = and ip and 
dip/dx to be continuous at p> = <p ma x- 

Our analysis throughout this paper focusses on the regime, where ion reflection is not 
present yet and provides the condition for ion reflection. An alternative way how to introduce 
dissipation and to break the soliton symmetry in this Mach number range, is taking into 
account Landau damping of electrons. The solitary electrostatic potential accelerates the 
electrons while propagating towards the peak and a considerable amount of electrons will 
have velocities around the resonance condition v e ~ v p h = c s . The interaction with ion 
acoustic waves leads to acceleration of particles with velocities v e < c s and deceleration of 
particles with velocities v e > c s . As the number of particles with velocities v e ^ c s is higher 
than the number of particles with velocities v e > c s , the waves will be damped. Landau 
damping of ions can become important if the ions are also trapped in the downstream 
potential. The consequence of this wave damping is a polynomial solution for the Sagdeev 



potential 



251 ] . which gives a shock-like solution for the electrostatic potential 24]. 
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VI. SUMMARY 



In this paper, electrostatic shock solutions have been identified starting from an initial 
system of two colliding plasma slabs, containing two populations of hot electrons and cold 
ions and a population of trapped electrons. From the calculation of the pressure terms, 
the Sagdeev potential has been derived from the initial conditions for arbitrary density 
and temperature ratios, and a range for the steady state Mach number was presented for 
the stage where ion reflection starts to become important. For the first time, relativistic 
electron temperatures have been considered and approximations for the highly relativistic 
case have been presented. By introducing the steady state bulk velocity of the downstream 
population, an actual dependence of the Mach number on the initial density and temperature 
ratios was gained, bridging the range between the minimum and maximum Mach numbers. 
The critical Mach number at which ion reflection appears, is achieved by increasing the 
electron temperature in the upstream plasma, which will lead to an increase of the sound 
speed c s and produce shock reflected ions with high energy. 

The solitary wave in the upstream is transformed into a shock solution if dissipative 
effects are taken into account that will break the symmetry. Such effects are for example 
the trapping of electrons and ions in the downstream, which can be theoretically described 
with the BGK approach, or the reflection of ions from the shock front, both described by 
a separate Sagdeev potential for the downstream. Also Landau damping can provide the 
necessary dissipation. The result is an oscillatory component in the downstream electrostatic 
potential. 
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